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Abstract.  The  first  documentation  of  Doppler  shifted 
auroral  Lyman  a  emission  resulting  from  incident  energetic 
protons  in  the  auroral  regions  has  been  made  using  nadir  VUV 
satellite  spectral  observations.  The  auroral  Lyman  a  emission 
from  high-velocity  protons  is  expected  to  show  a  red  shifted 
wavelength  displacement  based  on  ground-based  observations 
of  Balmer  lines.  VUV  spectra  (1 100-1900A)  taken  over  five 
sample  auroral  oval  crossings  by  a  nadir-viewing  satellite  in 
1978  consistently  show  the  Lyman  a  emission  displaced  to¬ 
ward  longer  wavelengths  with  a  larger  line  width.  The  inten¬ 
sity  peaks  were  shifted  up  to  4  A  when  the  geocoronal  Lyman 
a  emission  profile  was  subtracted  from  the  Lyman  a  emission 
profile  observed  over  the  auroral  regions.  The  optical  observa¬ 
tions  infer  the  auroral  proton  precipitation  with  average  eper-  / 
gies  of  34  keV  and  an  energy  flux  of  0. 1  erg/cm  7s ’ 4r  ' 
when  interpreted  according  to  available  model  calculations. 
These  values  agree  reasonably  well  with  the  average  values 
for  the  characteristics  of  nightside  incident  auroral  protons 
based  on  previous  statistical  satellite  panicle  precipitation  ob¬ 
servations. 

Introduction 

Doppler  shifted  hydrogen  Balmer  series  emission  in  the 
auroral  regions  was  first  observed  by  a  ground-based  spec¬ 
trometer  and  used  to  infer  the  energy  of  precipitating  protons 
[Meinel,  1951).  Subsequent  calculations  and  discussions  have 
been  presented  by  Chamberlain  [1961 )  and  Eather  [1967],  Fig¬ 
ure  1  shows  a  Doppler  shifted  and  broadened  Balmer  0  line 
observed  in  a  proton  aurora  compared  to  the  line  observed 
in  the  laboratory  (Zwick  and  Shepherd,  1963] .  Auroral  pro¬ 
tons,  typically  10-40  keV,  precipitate  into  the  atmosphere  and 
lose  energy  through  charge-exchange  and  charge-stripping  cy¬ 
cles  with  atmospheric  constituents.  Relatively  recent  detailed 
computations  of  the  profile  of  the  hydrogen  spectral  line  emis¬ 
sion  were  performed  for  ground-based  observations  in  the 
directions  of  the  magnetic  zenith  and  horizon  [Ponomarev, 
1976],  The  profile  showed  a  good  match  to  the  observed 
Balmer  line  spectra  [Galperin  et  al.,  1976]. 

The  Lyman  a  emission  line  <  1216  A )  from  proton  precipita¬ 
tion  can  be  observed  only  from  space.  In  this  study  we  present 


the  first  observations  of  the  Doppler  shifted  profile  of  the 
auroral  Lyman  a  line  deduced  from  a  polar  orbiting  satellite 
and  infer  the  average  energy  and  energy  flux  of  the  auroral 
protons  by  comparing  the  observed  Doppler  shift  and  broad¬ 
ening  with  the  available  model  calculations.  The  Lyman  a 
Doppler  shift  of  10-eV  hydrogen  atoms  on  Jupiter  was  re¬ 
cently  reported  in  the  American  Astronomical  Society  meet¬ 
ing,  Division  of  Planetary  Sciences  [Clarke  et  al.,  1988],  as 
pointed  out  by  the  reviewer  of  this  paper. 

Observations 

The  auroral  spectral  data  were  obtained  by  an  AFGL  UV 
background  experiment  flown  on  the  S3-4  satellite  in  1978 
[Huffman  et  al.,  1980].  On  board  the  satellite,  which  was  in  4 
a  low-altitude  (180-270  km)  polar  orbit  near  the  noon-mid- 
night  meridian  plane,  was  the  nadir-viewing  1/4-m,  f/5  Ebert- 
Fastie  spectrometer.  Fifteen  spectra  with  8  A  resolution  were 
selected  for  analysis  during  periods  of  strong  Lyman  emis¬ 
sion.  All  the  spectra  were  taken  at  a  satellite  altitude  of  about 
260  km.  Table  1  lists  the  geometrical  and  geophysical  infor¬ 
mation  along  with  the  spectral  characteristics  from  five  auroral 
oval  crossings.  Figure  2a  is  an  example  of  the  nadir  Lyman 
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Fig.  1  Examples  of  hydrogen  Balmer  0  line  profiles  [from 
Zwick  and  Shepherd,  1963);  (a)  a  hydrogen  lamp  source,  (b) 
a  magnetic  zenith  profile  of  a  proton  aurora. 
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TABLE  1.  Summary  of  Lor  spectra. 
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a  line  intensity  as  observed  along  the  entire  orbit  on  May  1 , 
1978.  The  geocoronal  Lyman  a  line  intensity  shows  a  strong 
correlation  with  the  solar  zenith  angle.  Figure  2b  shows  two 
spectra,  one  from  the  proton  auroral  region  and  the  other 
from  the  midlatitude  region.  The  latter  is  the  average  of  22 
nighttime  geocoronal  spectra  to  reduce  noise  and  has  the  8  A 
half  width  full  maximum  (HWFM)  line  profile.  The  former 
has  a  broader  10A  HWFM  spectral  profile  with  a  peak  red 
shifted  by  2A,  indicating  that  some  emission  is  from  high¬ 
speed  hydrogen  moving  away  from  the  satellite  toward  the 
earth. 

Two  steps  were  needed  in  order  to  deduce  the  auroral  Ly- 
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Fig.  2  (a)  Observed  Lyman  a  line  intensity  and  solar  zenith 
angle  of  the  S3-4  throughout  an  orbit,  (b)  Comparison  of 
the  Lyman  a  profiles  taken  from  auroral  and  midlatitude 
regions  (the  corresponding  observation  times  were  indicated 
in  Figure  2a). 


man  a:  first,  the  geocoronal  background  intensity  was  esti¬ 
mated  by  interpolation  of  the  smoothed  geocoronal  intensities 
outside  the  auroral  region;  second,  the  geocoronal  Lyman  a 
profile  with  the  estimated  background  intensity  was  subtracted 
from  the  observed  Lyman  a  profile  over  the  proton  auroral 
region  as  shown  in  Figures  3a-f.  The  peak  shifts  and  the 
HWFM  width  of  the  remaining  auroral  spectra  are  shown 
in  Table  1;  the  auroral  Lyman  a  spectrum  is  broader  than 
that  of  the  geocorona  by  -2A. 

Discussion  and  Conclusions 

We  used  the  hydrogen  line  spectral  profiles  of  auroral  pro¬ 
tons  computed  by  a  Monte  Carlo  method  as  our  model 
[Ponomarev,  1976].  The  model  takes  into  account  the  charge- 
exchange  and  proton-beam  stripping  processes  in  the  three 
dimensional  model  atmosphere  and  uses  the  cross  sections 
of  McNeal  and  Birely  [1973].  Figure  4  shows  the  profile  cal¬ 
culated  for  an  isotropic  pitch  angle  distribution  as  a  func¬ 
tion  of  the  Doppler  velocity  in  two  viewing  directions.  The 
centroid  of  the  profile  shifts  to  a  more  negative  Doppler  ve¬ 
locity  as  seen  from  the  ground  with  increasing  energy  of  the 
monoenergetic  beam.  In  order  to  compare  this  model  with 
our  space-based  observations,  the  profiles  of  Figures  4b  and 
c  were  transformed  for  the  1216  A  line  and  also  smeared  for 
8  A  spectral  resolution.  Figure  5a  shows  the  increase  of  both 
the  Doppler  shift  and  the  broadening  with  the  increasing  ener¬ 
gy  of  the  protons  as  viewed  downward  from  the  satellite.  In 
order  to  reduce  the  noise  level,  we  averaged  six  spectra  from 
the  same  auroral  oval  pass  on  July  5  and  subtracted  the  esti¬ 
mated  geocoronal  Lyman  a  profile  (Figure  5b).  Figure  5c  il¬ 
lustrates  the  remaining  (i.e. ,  auroral)  spectrum  and  two  model 
spectra  calculated  for  two  proton  energies.  The  general  shape 
of  the  auroral  profile  matches  the  calculated  30-keV  profile. 
The  observed  profile  is  broader  than  the  30-keV  model  profile, 
as  expected,  because  the  energy  spectrum  of  the  auroral  pro¬ 
tons  is  not  a  monoenergetic  beam  and  the  deduced  auroral 
spectrum  is  from  an  average  of  six  spectra.  Taking  these  facts 
into  consideration,  the  auroral  spectrum  matches  the  30-keV 
model  calculated  profile  very  well. 

According  to  Ponomarev’s  calculation  for  isotropic  proton 
precipitation,  the  spectral  peak  shift,  A\,  can  be  related  to 
the  proton  energy,  E0,  by  the  equation; 
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Fig.  3.  Examples  of  Lyman  a  profiles  from  an  auroral  region.  The  estimated 
geocoronal  Lyman  a  profile  (dotted  line)  was  subtracted  from  the  observed 
profile  to  reveal  the  Lyman  a  profile  from  the  proton  precipitation 
(designated  (1)  —  (2)). 
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E„  =  12.321  x  AX  -  6.291. 

Using  this  equation,  the  mean  energy  of  auroral  protons  can 
be  estimated  (Table  1).  The  values  are,  in  general,  higher  than 
those  from  statistical  study  of  ion  precipitation  measurements 
reported  by  Hardy  et  al.  [1988],  which  are  about  20  keV  in 
the  midnight  region  for  Kp  >  3.  However,  the  DMSP  mea- 


Fig.  4.  Hydrogen  line  profile  in  the  proton  aurora  as  a  func¬ 
tion  of  the  Doppler  velocity  computed  by  the  Monte  Carlo 
technique.  The  geomagnetic  zenith  (7.M)  and  horizon  (HM) 
profiles  are  calculated  for  an  isotropic  distribution  of  the  pitch 
angle  for  protons  of  monoenergetic  energies  of  1,  5,  and  30 
keV  (from  Calperin  ct  al.  1976j. 


surement  was  limited  to  energies  from  30  eV  to  30  keV  and 
some  spectra  were  extrapolated  to  100  keV  to  include  the  high- 
energy  tail  of  the  ion  distribution.  This  extrapolation  increased 
the  average  energy  and  the  energy  flux  by  less  than  30%  and 
50%,  respectively.  The  Lyman  a  spectra  used  here  were 
among  the  most  intense  proton  emissions  in  our  satellite  data 
set,  corresponding  to  very  active  periods  (up  to  Kp  =  6- ). 
A  higher  mean  energy  than  the  statistical  value  is  expected, 
because  there  is  no  temporal  or  spatial  averaging,  which  would 
smooth  down  the  high  values.  Therefore,  the  agreement  be¬ 
tween  our  inferred  and  averaged  measured  proton  energy  is 
reasonable. 

For  the  energy  flux  estimation,  we  used  the  Lyman  line 
emission  intensity  calculated  by  Edger  et  al.  (1973)  with  the 
emission  cross  section  increased  by  22%  as  suggested  by  Van 
Zyl  and  Newmann  [1988],  The  estimated  energy  fluxes  are 
derived  using  the  proton  energy  inferred  from  the  observed 
profile  (the  last  column  of  Table  1).  The  average  vaJue  esti¬ 
mated  from  the  observation  is  -  0.3  erg  cm  :  s'1.  This  val¬ 
ue  would  be  divided  by  3  for  the  integrated  solid  angle  for 
comparison  with  the  statistical  particle  precipitation  measure¬ 
ments  of  0.1  erg  cm  2s  !sr  1  reported  by  Hardy  et  al. 
(1988}.  Thus,  the  observed  auroral  Lyman  a  profile  and  emis¬ 
sion  intensity  are  consistent  with  the  existing  model  calcula¬ 
tions  and  ion  precipitation  measurements.  These  observations 
suggest  that  sensitive  high-resolution  spectral  measurements 
from  a  satellite  can  be  used  for  global  imaging  of  the  charac¬ 
teristics  of  the  proton  precipitation.  However,  specific  cal¬ 
culations  of  the  resulting  Lyman  a  emission  profile  expected 
for  different  nadir  viewing  angles,  including  multiple  scat¬ 
tering  effects  on  the  emission  profile,  would  be  required. 


□  □ 


146 


Ishimoto  et  al.:  Doppler  Shift  of  Auroral  Lyman  a 


i.o- 

0.8 

|  0.6 

t  0.4 
c 
a> 
c 

“  0.2- 


(•> 


En  =  5  keV 


6  ,->  =  30  keV 


1205  1210 


.  1  . 1  .  .  _ 

1216  1220 
Wavelength  (A) 


1225  12*3  cT  ?235 


200 


(b) 


160 


< 

i  120 


80 

40 


Observed 


Geocorona 


0  '  .Tr; 

1205 


'  1210 .  1 216  1220  i  225~  ~V22Q 

Wavelength  (A) 


1235 


1.0 - 

.  (c) 

0.8  E0  ■  5  keV 

> 

I  0.6 


Eg  *  30  keV 

Aurora 


0  - — 

1205 


j . .{.. . J . i._ 

1210  1216  1220  1225 

Wavelength  (A) 


1230  1235 


Fig.  5.  Comparison  of  the  auroral  Lyman  a  profiles:  observed 
and  calculated,  (a)  Lyman  a  spectral  intensity  profile  in  the 
proton  aurora  adapted  from  model  calculations  of  Ponomarev 
[1976].  (Figures  4b  and  4c  are  translated  from  the  Doppler 
velocity  to  the  wavelength.) 

(b)  Auroral  Lyman  a  spectral  intensity  profile  observed  on 
June  5,  1978.  OBSERVED:  an  average  of  six  consecutive 
auroral  spectra.  GEOCORonia-  he  estimated  intensity  with 
a  profile  of  an  average  o  J  rn  cilatitude  spectra  from  the 
same  orbit.  AURORA:  OP  ,iD- GEOCORONA. 

(c)  Comparison  of  observed  at  :al  Lyman  a  spectral  rela¬ 
tive  intensity  profile  (AURORA  profile  in  (b))  and  model  cal¬ 
culation  (in  (a)). 
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